A Ball in a Groove 
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We study the static equilibrium of an elastic sphere held in a rigid groove by gravity and frictional 
contacts, as determined by contact mechanics. As a function of the opening angle of the groove and 
the tilt of the groove with respect to the vertical, we identify two regimes of static equilibrium for 
the ball. In the first of these, at large opening angle or low tilt, the ball rolls at both contacts as it 
is loaded. This is an analog of the "elastic" regime in the mechanics of granular media. At smaller 
opening angles or larger tilts, the ball rolls at one contact and slides at the other as it is loaded, 
analogously with the "plastic" regime in the mechanics of granular media. In the elastic regime, the 
stress indeterminacy is resolved by the underlying kinetics of the ball response to loading. 

46.55.+d, 45.70.Cc, 46.25.-y 

The behavior of granular materials is one of the most 
difficult subjects in solid mechanics Q. It is elemen- 
tary to state the problem-given a packing of elastic bod- 
ies, whose contacts with one another obey the classical 
Amonton's laws of friction, what is the nature of the 
stress transmission in the system in the limit where the 
forces between the bodies are much smaller than the 
appropriately dimensioned elastic moduli of the bodies? 
This problem underlics much of soil mechanics, sedimen- 
tology, and powder technology, but has thus far resisted 
a priori solution, although much is known at an empirical 
level about the behavior of such systems @. 

The difficulty of this problem ariscs from its non- 
linearity, which appears in two different guises. In the 
first place, the forces between elastic bodies in contact 
are non-linear functions of their displacements. In the 
second case, sufficiently strong shear forces can cause a 
contact to fail, which is also a source of non-linearity. 

Soil mechanicists often focus on the second source of 
non-linearity, but presume that at a macroscopic level 
the first source can be hnessed. The result is a class of 
"elasto-plastic" models that focus on the interplay be- 
tween an assumed linear elastic region of soil response 
and the plasticity connected with shear-induced failurc 
Whilc this is a plausible and effective approach, it 
is disturbing to realize that there is no fundamental rea- 
son to expect linear elastic behavior in a granular system 
at all, given the non-linearity of contact mechanics. This 
uncertainty in the fundamental underpinnings of gran- 
ular statics is sometimes termed the "stress indetermi- 
nacy" problem, as the stress tensor cannot be determined 
from the equations of force balance alone; these must be 
supplemented with some constitutive relations 

This uncertainty has motivated statistical approaches 
to the problem, notably the "q -model" and its variants, 
which focus on the fluctuations of forces within a granular 
packing §| . Unfortunately, these models are not based on 
a realistic picture of the contact forces between particles. 

For this reason, we here address a simple model sys- 
tem, which nevertheless exhibits stress indeterminacy: a 



FIG. 1. A ball in a groove is subject to five forces-the gravi- 
tational force, which acts on its center of mass, and tangential 
and normal forces at each contact. The two rigid planes are 

inclined at angles 6, l . to the horizontal. 
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ball placed in a groove with perfectly rigid walls (Figure 
|ÏJ) under gravitational loading. Since the contact me- 
chanics of this system are tractable, we are able to solve 
its stress indeterminacy problem. We find that it is pos- 
sible to specify an "elastic" and a "plastic" regime in the 
behavior of this ball, although the elastic regime is not 
closcly rclatcd to ordinary continuüm elasticity. In the 
elastic regime, the stress indeterminacy is lifted by the 
underlying kinetics of the ball as it is loaded-since the 
number of kinematic parameters in this case is identical 
to the number of constraints on the ball, we term this 
the isokinetic regime, in contrast to the isostatic plastic 
regime Q, in which the number of independent forces 
equals the number of constraints. 

Suppose that, as in Fig. |ï|, we have an elastic sphere 
of radius R, Young's modulus E, Poisson ratio v, and 
mass M in contact with a rigid groove. We consider only 
motions in the plane perpendicular to the groove axis 
(out of the page in Fig. [ï]), thus reducing our problem to 
an essentially two dimensional one (note however, that 
we still consider an elastic sphere-the contact mechanics 
of a disk are essentially different.) Or, and Or denote the 
angle between the left and right walls of the groove with 
respect to the horizontal, respectively. Without loss of 
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generality we choose 9l > 9 R . It is convenient to define 
7 = + 0r) and S = ~ &r), so that n — 2j is the 
groove opening angle, and 6 gives the tilt of the groove 
with respect to the vertical. Our definition of a groove 
satisfics < 6 < 7 < § . The normal and tangential 
forces exerted by the wall upon the ball are Nl, Tl on 
the left, and N R , T R on the right, with the tangential 
forces defined so that a force exerted upwards is positivo. 
The ball is also subject to the gravitational force Mg. 
The static equilibrium of the ball requires that 

T L sin L + N L cos 9 L + T R sin 9 R + N R cos 9r = Mg, (1) 
-T L cos 9 L + N L sin 6 L + T R cos 9 R - N R sin 9 R = 0, (2) 

Tl = T R , (3) 

corresponding to the requirements that the net horizon- 
tal and vertical forces, and the net torque, be zero. In 
addition, Amonton's law requires that 

r (S) <*tf (S) , (4) 

where k is the coëfficiënt of friction between the ball and 
the wall. Since we have three constraints for the four 
forces Tl, Nl, T r , and Nr, even subject to the Amon- 
ton's law inequalities we will have in general a manifold 
of solutions. In order to determine which solution on this 
manifold is actually chosen we must examine in more de- 
tail the mechanics of the ball itself. 

The problem of an elastic sphere compressed against a 
rigid flat surface was given its classical solution by Hertz 
. If the sphere (of radius R) is loaded by a force normal 
to the wall N, then its center will displace towards the 
wall by a distance u, forming thereby a contact circle of 
radius a, with 

a = (GNR)i, (5) 

where G = 3(1 - v 2 )/(AE). 

The situation in the presence of tangential forces is 
considerably more complicated 0. If an already estab- 
lished normal contact is sheared at a tangential force that 
is insufficiënt to cause failure, a region of microslip adja- 
cent to the outer perimeter of the contact forms, moving 
inwards as the tangential force increases. In this region, 
the two surfaces slip with respect to one another. In the 
interior of the contact, the two surfaces remain stuck to- 
gether as the tangential force is applied. Eventually, at 
failure, there is no stick region left. 

The details of this process were explored in a classic pa- 
per by Mindlin and Deresiewicz ||. They observed that 
the analysis of the contact structure at finite tangential 
and normal force is only tractable in the limit of low Pois- 
son ratio v. In this limit the local normal and tangential 
strains at the contact decouple from one another, so that 



the normal pressures at the contact are always those of 
a Hertzian contact with the same overall normal force. 
Thus the additional tangential strains do not alter the 
solution of the part of the problem corresponding to the 
normal force. 

The problem is essentially hysteretic, in that the 
strains corresponding to a given tangential and normal 
force are determined by the specific history of contact 
loading leading to these forces. A general increment of 
force in such a loading history can be described by the 
change in tangential force ÖT and the change in normal 
force 6N. A special class of such increments do not in- 
volve microslip over some area of the contact surface, 
those with || 

SN > 0, (7) 
\ST\ < kSN. (8) 

This latter result suggests that we specify a particu- 
larly simple loading history of a contact, in which T/N 
is maintained constant, with an absolute value < k, while 
one force, say N, is monotonically (and adiabatically) in- 
creased as a function of time. This type of loading history 
we term proportional loading. 

For proportional loading, we can also specify the com- 
pliance of the contact, or the values of u n , Ut, which are 
respectively the normal and tangential displacements of 
the center of the ball corresponding to fixed forces N and 
T. From the Mindlin-Deresiewicz results, it is straight- 
forward to show that for small Poisson ratio, 



with u n given by the Hertzian result, Eq. (|^). 

There are other types of motions possible for a contact. 
One is torsional motion, corresponding to the rotation of 
the ball about an axis normal to the plane of the contact. 
Torsional motion leads to couples that vanish as the con- 
tact size goes to zero, so we will disregard this type of 
motion. Rotation of the partiele about an axis parallel 
to the plane of the contact, on the other hand, leads to a 
rolling motion, in which the center of the contact moves 
parallel to the underlying plane. 

Unfortunately, the strains corresponding to rolling mo- 
tion have not been successfully computed analytically, 
even for small Poisson ratio, although a number of ap- 
proximate solutions have been proposed. For rolling at 
fixed normal force, there is always a region of microslip 
at the trailing edge of the contact 0. However, it does 
not necessarily follow that a contact that rolls while be- 
ing proportionally loaded experiences microslip, since the 
normal force will be increasing in this case. In any case, 
proportional loading will lead to a unique relation be- 
tween Ut/u n and T/N, where Ut is defined as the trans- 
verse displacement of the ball center of mass with respect 
to the contact, not with respect to its initial position. For 
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our purposes, it will be sufficiënt to use Eq. (j^j as an ap- 
proximate form, since the qualitative properties of the 
solution are relatively independent of the exact form of 
the compliance. 

In order to insure that the contacts are proportion- 
ally loaded, we consider the following experiment. First 
the ball is placed into the groove at zero gravity so that 
it barely touches both walls of the groove. Then the 
gravitational constant is slowly increased from zero to its 
physical value. Since the external force during this load- 
ing increases and does not change direction, we expect 
T oc N at each contact with N monotonically increasing. 

Let us suppose that the center of mass of the ball moves 
a distance Z along a path that makes an angle Bl with 
the left-hand wall and /3r with the right-hand wall (Note 
that Bl + Pr = n — 27.) During the loading, the ball 
may also rotate as a rigid body by an angle <f>, with </> > 
denoting counterclockwise rolling. Finally, the position 
of the contacts with the wall will in general move by 
M c .( £ ) on the left and right-hand walls respectively; u c > 
corresponds to contact motion down the wall. 

Now we must specify the normal and tangential dis- 
placements of the center of mass of the ball with respect 
to the contact positions on the two walls. The normal 
displacements u are relatively simple to write: 



(10) 



The tangential displacements are a little bit more com- 
plicated, because of rolling, and because the left and 
right-hand contacts may have slipped by l ^ during 
the loading, respectively. Then, u c -l = —Zv + ó,l and 
u c -r = Zv + d,R, where v = R<p/Z. Thus, the tangential 
displacements u t .^L^ are given by 



1. , l = Ztcos 8, ±v) — d, l-,. 



We are now using the fivc kinematic variables Z , Bl, v, 
and d, l . . The proportional loading assumption insures 
that whether the contacts are loaded with or without 
rolling, the forces still obey 



T 



(12) 



where Amonton's law insures that this is true for sliding 
contacts as well. 

Rolling solutions We first consider solutions in which 
both contacts roll without sliding: One of the contact 
points retrogrades up its contact plane, while the other 
rolls down its plane. Thus we have Eqs. (10) for the 
normal displacements, while Eqs. ([Dl) become 



l t;( R ) = Z{cos3 {Ir) ±v), 



In principle, we have a uniquely soluble system, since 
we have the three variables Z, Bl and v describing the 



kinematics, with three conditions of static equilibrium. 
However, we must also impose Eqs. (|j) constraining the 
values of T/N. Thus there will be some restricted region 
of the 7 — 5 plane in which solution of the Eqs. (|ï[]|) 
is possible. Substituting Eqs. ( JTÖj ) and ( |ï"3| ) into Eqs. 
(|ïl||), we obtain 



3 

sin Br \ 2 sin ^ (cos 8l + v) (cos 6r — cos 



sin/3 L 



sin Br 



sin Bl sin Or 



and 



1 

'sin/JfiA 2 cos Bl 



- v 



v sin B L ) cos Br - 



v 



1) 



(14) 



(15) 



It is straightforward to solve these equations for d L 
(and thus Br) and for v. For 0l = Or (or S = 0), Bl = 
Br = 77 — 7 and v — 0, i.e., there is no rolling and both 
contacts are stuck. Furthermore, 



Uf.L 



Uf,R 



= cot0L = tan7 = = — ^ < k. (16) 



Un-L Un;R. 



N L N r 



Hence, the solution for 5 = can only exist for 7 < 
arctan(/c), because for 7 > arctan(fc), the values of T/N 
corresponding to this solution exceed k. Thus, emanat- 
ing from the point (7,(5) = (arctan(fc), 0) is a bound- 
ary beyond which doublé rolling is impossible, because 
it would require illegal values of T/N > k. This bound 
is exceeded first on the steeper of the two walls. This 
boundary, shown as a solid line in Fig. ^| rejoins the 
7-axis at 7 = 0, asymptotically merging with the line 
5 = 7. 

There is an additional solution for which v = (no 



( n ) rolling), with B L ^ Pi- 



that exist s for j 



< 



< 



arctan(v / 2). This additional solution always has at least 
one of the ratios T/N > y/2, so it can only appear for 
k > v2> and constitutes a boundary between a regime 
in which the ball rolls clockwise (i.e., down the steeper 
surface) for smaller 7, and a regime in which the ball rolls 
counter-clockwise (i.e., up the steeper surface) for larger 
7. The values of (7, S) for which this solution occurs is 
shown for k = 10 as a dot-dashed line in Fig. 

Rolling-sliding solutions Beyond the boundary at 
which the steeper wall has T = kN, the only solutions 
possible are "rolling-sliding" solutions, in which one con- 
tact rolls without sliding while the other slides. Assuming 
that the left-hand wall slides, we have 



N L 

(13) The equations of static equilibrium then imply 



T _ smÖR 
N R coa 0r - cos 0l + \ sinöi ' 



(17) 



(18) 



3 




90 



60 



30 



30 60 90 

Y (degrees) 

FIG. 2. Boundaries between different regimes of behavior 
for coëfficiënt of friction k — 1 (thin lines) and k = f0 (thick 
lines), as a function of 7 = + 9r) and 8 = — 0r). 
Solid lines mark the boundary between the elastic dou- 
blé-rolling regime below and the plastic rolling-sliding regime 
above. The dashed and dot-dashed lines mark the boundaries 
between clockwise (low-7) and counter-clockwise (high-7) 
rolling. 



where fi depends only upon the geometry and upon the 
coëfficiënt of friction k. Note that for 9r < 9l, /•* < k, 
while for 9r > 9l, fi > k. Thus, sliding always occurs on 
the steeper wall, which is indeed the left-hand wall for 
6 > 0. This is the simplest continuation of the double- 
rolling solution, which failed at the boundary through 
incipient sliding at the steeper wall. We have 



Nl 
N r 



smfi L 
sin/3 fi 



(19) 



which can be uniquely solved for Pl and Pr. In addition, 
the criterion that the ball rolls at the right-hand contact 
implies that 



T 

Nr~ 



/' 



COS, [3r 



sin Pr 



(20) 



which allows us to solve for v(J3r). 

The physical constraint > sets the boundary 
on the region with a rolling-sliding solution, which, for 
7 < arctan(fc), coincides with the upper boundary of 
the doublé rolling solution discussed in the previous sec- 
tion. For 5 = we have v = and T/Nr — k, and 
this solution collapses onto a doublé sliding solution for 
7 > arctan(fc). In addition, there is a line of Solutions 
with v — 0, corresponding to one contact sliding and 
the other being stuck. For k < \/2 this line emanates 
from (7,0) = (arctan(fc), 0) and ends at (7, S) — 
. For k > -\/2, this line emanates from the point on the 
boundary of the doublé rolling region for which v = 0, 
and ends at (7,(5) = (f-, j) (dashed lines in Fig. . As 



for the doublé rolling v = line, this line separates the 
low-7 clockwise rolling from the high-7 counterclockwise 
rolling. 

Much discussion of the statics of granular media has 
cmphasized the two limits of plastic and linear elastic be- 
havior. Clearly the rolling-sliding regime of the ball in 
the groove is a good analog to plastic behavior-in this 
regime, the problem is isostatic, meaning that the num- 
ber of independent forces is identical to the number of 
constraints. In the doublé rolling regime, the state of the 
ball is specified by its displacements, as in conventional 
elasticity theory. However, minimization of the strain en- 
ergies corresponding to the contact forces plays no role 
in resolving the stress indeterminacy. 

It is instructive to consider the relationship between 
the forces and the underlying degrees of freedom Z, Pl , 
v. and d. 



*(£)• 



The four forces T, l , and N, l ■, are functions 

Kr) Kr) 



of these five underlying variables. In the double-rolling 



regime, 



L ( L n) 



0, and there are only three active vari- 



ables, which are then uniquely fixed by the equations of 
static equilibrium. We call this an isokinetic configura- 
tion, because the number of kinetic degrees of freedom 
equals the number of constraints enforcing equilibrium. 

On the other hand, as one passes into the isostatic 
rolling-sliding regime, one acquires a new degree of free- 
dom d/,, so that there are now four underlying active 
degrees of freedom. However, one of these degrees of 
freedom is used to fix Tl = kN^, so in practice there are 
three degrees of freedom remaining to satisfy the equa- 
tions of static equilibrium. 

We are grateful to P.M. Chaikin and D. Levine for 
stimulating discussions. 
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